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A TROPICAL ATMOSPHERE MODEL WITH MOISTURE: 
GLOBAL WELL-POSEDNESS AND RELAXATION LIMIT 


JINKAI LI AND EDRISS S. TITI 

Abstract. In this paper, we consider a nonlinear interaction system between the 
barotropic mode and the first baroclinic mode of the tropical atmosphere with 
moisture; that was derived in [Frierson, D. M. W.; Majda, A. J.; Pauluis, O. M.: 
Dynamics of precipitation fronts in the tropical atmosphere: a novel relaxation limit, 
Commum. Math. Sci., 2 (2004), 591-626.] We establish the global existence and 
uniqueness of strong solutions to this system, with initial data in , for each fixed 
convective adjustment relaxation time parameter e > 0. Moreover, if the initial data 
enjoy slightly more regularity than H^, then the unique strong solution depends 
continuously on the initial data. Furthermore, by establishing several appropriate 
e-independent estimates, we prove that the system converges to a limiting system, 
as the relaxation time parameter e tends to zero, with convergence rate of the 
order 0{\/e). Moreover, the limiting system has a unique global strong solution, 
for any initial data in and such unique strong solution depends continuously 
on the initial data if the the initial data posses slightly more regularity than H^. 
Notably, this solves the VISCOUS version of an open problem proposed in the 
above mentioned paper of Frierson, Majda and Pauluis. 


1. Introduction 


1.1. The primitive equations for planetary atmospheric dynamics. In the 

context of large-scale atmosphere, the ratio of the vertical scale to the horizontal 
scale is very small, which, by scale analysis, see, e.g., jss, 42 \, leads to the hydrostatic 
approximation in the vertical momentum equation. This small aspect ratio limit can 
be rigorously justified, see [I, 29|. Taking into account the Boussinesq approximation 


and the hydrostatic approximation to the Navier-Stokes equations, one obtains the 
primitive equations, which model the large-scale atmospheric dynamics. 

The primitive equations read (see, e.g., 0 0 S m 0 Hi, Ell) 

dtV + (V ■ V^V + Wd,V - pAV + V,,<I> = 0, 

= f ^ 

dt© + V ■ V,.© + WdS + 

Wh-Y + d,W = Q, 


= S, 

9 


0 , 


( 1 . 1 ) 
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where the unknowns V = (Vi,V 2 )^, W, $ and 0 are the horizontal velocity held, 
vertical velocity, pressure and potential temperature, respectively, while the positive 
constant /i is the viscosity coefficient. The total potential temperature is given by 

0total(^, 


where 6q is a positive reference constant temperature and 6 dehnes the vertical prohle 
background stratihcation, satisfying = {g/9o)dz9 > 0, where N is the Brunt- 
Vaisala buoyancy frequency. Here we use V/j = (d^ydy) to denote the horizontal 
gradient and V"*" = (—V 2 , Hi)^. 

During the last two decades, a lot of efforts have been done on the mathematical 
studies of the primitive equations. Up to now, it has been known that the primitive 


equations, with full viscosity and full diffusivity, have global weak solutions (but the 


and have a unique global strong solution, 
sTIgI. Il2l. 2^ for some recent developments towards 


uniqueness is still unclear), see 
see [ill. 24, [ 2 ^, and also see 
the direction of partial dissipation cases. Moreover, the recent works [7|-|9| show that 
the horizontal viscosity turns out to be more crucial than the vertical one for the 
global well-posedness, because the results there show that the vertical viscosity is not 
required for the global well-posedness of strong solutions to the primitive equations. 
Notably, the invicid primitive equations may develop hnite time singularities, see 
441. Combining the results of [3^ and those of 01 e can conclude that 


the horizontal viscosity is necessary for the global well-posedness of the primitive 
equations, and if ignoring the temperature effect, the horizontal viscosity is also 
sufficient for the global well-posedness. 


1.2. The barotropic and the first baroclinic modes interaction system. In 

the tropics, the wind in the lower troposphere is of equal magnitude but with opposite 
sign to that in the upper troposphere, in other words, the primary effect is captured 
in the hrst baroclinic mode. However, for the study of the tropical-extratropical 
interactions, where the transport of momentum between the barotropic and baroclinic 
modes plays an important role, it is necessary to retain both the barotropic and 
baroclinic modes of the velocity. 

Consider the primitive equations fll.ip in the layer x {0,H), for a positive 
constant H. Since we consider the tropical atmosphere and take into consideration 
the tropical-extratropical interactions, we can impose an ansatz of the form 

I ix.v.z.t) =1 

pj y ,U, , 


^ I = {'^) {x, y,t)+ [^ ) {x, y, t)V2cos{7iz/H) 


and 


{x, y,z,t)= ( ^ ) {x, y, t)V2sm{nz/H), 


{x,y,z,t) = 

which carry the barotropic and hrst baroclinic modes of the unknowns. 

By performing the Galerkin projection of the primitive equations in the vertical 
direction onto the barotropic mode and the hrst baroclinic mode, one derives the 
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following dimensionless interaction, between the barotropic mode and the first baro- 
clinic mode, system for the tropical atmosphere (see 
the details): 


and also 15, 19, 35, 41 for 


{ dtu + {u ■ V)m — Am + Vp + V ■ (n ® n) = 0, 
V • M = 0, 

dtv + {u ■ V)m — Am + {v ■ V)m = V6*, 

dte + u-Ve-V-v = Se, 


where u = (mi,M 2 ) is the barotropic velocity, and v = (mi,M 2 ), p and 0, respectively, 
are the first baroclinic modes of the velocity, pressure and the temperature. The 
system is now defined on M^, and the operators V and A are therefore those for the 
variables x and y. 


1.3. The moisture equation. An important ingredient of the tropical atmospheric 
circulation is the water vapour. Water vapour is the most abundant greenhouse gas 
in the atmosphere, and it is responsible for amplifying the long-term warming or 
cooling cycles. Therefore, one should also consider the coupling with an equation 
modeling moisture in the atmosphere. 

Following [l^ , we couple system fll.2p with the following large-scale moisture equa¬ 
tion 

dtq + u-Vq + QV -v =—P, (1.3) 

where Q is the prescribed gross moisture stratification. The precipitation rate P is 


parameterized, according to [15|, l20|, 1^, l41 


as 


P=-{q 


a9 — q)~ 


(1.4) 


where /■•■ = max{/, 0} denotes the positive part of /, £ is a convective adjustment 
time scale parameter, and a and q are constants, with g > 0. 

In order to close system fll.2l) - fll.3p . one still needs to parameterize the source term 
Se in the temperature equation. Generally, the temperature source Se combines three 
kinds of effects: the radiative cooling, the sensible heat flux and the precipitation P. 
For simplicity. 


and as in 
source term, i.e., we set 

with P given by fll.4p . 


36| . we only consider in this paper the precipitation 


Se = P, 


As in 15, 3^, by introducing the equivalent temperature Tg and the equivalent 


as 


moisture q^ 

Te = q + 9, 

system fll.2p - fll.3ll can be rewritten as 


= q — a9 


qe = q 


dtU + {u ■ V)m — Am -|- Vp -I- V • (m 0 m) = 0, 
V ■ M = 0, 


(1.5) 

( 1 . 6 ) 
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(1.7) 

( 1 . 8 ) 

(1.9) 


dtTf, + u ■ VTe — (1 — Q) V • V = 0, 



in X (0, oo), where the constants a and Q are required to satisfy (see 0 ) 


0<Q<1, a + Q > 0. 


( 1 . 10 ) 


1.4. Main results. We will work in the framework of strong solutions, which are 
dehned below. 

Definition 1.1. Given a positive time T and the initial data (mq, vo, 7"e,o, Q'e,o)- A 
function {u, v, T^, qf) is called a strong solution to system U.5\) - n~^) . on x (0, T), 
with initial data (mq, Uoi ^e, 0 i ?e,o)i enjoys the following regularities 


{u,v) e C{[Q,r\-,H\M}))nL\Q,T- 

{dtU,dtV,dtT,,dtq,) e L2(0,T;L2(M2))^ 

(Te, qe) e ^([0, 71; L2(M2)) n L-(0, T; 


and satisfies equations / Ij.5l) -/ f773) . a.e. on x (0,T), and has the initial value 


(m, n, T^e, (Je) |t=0 (iIQ WO) ^e,0) *?e,o) • 


Definition 1.2. A function {u,v,Tg,qe) is called a global strong solution to system 
/ [i.5)) -/ f77P)) . if it is a strong solution to system on x (0,T), for any 

positive time T■ 

Throughout this paper, for positive integer k and positive q G [1, cxd], we use 
and to denote the standard Lebesgue and Sobolev spaces, respectively, and 

when g = 2, we use instead of 1T*^’^(M^). For simplicity, we usually use \\f\\q 

to denote the ||/|| 2 , 9 (r 2 ). 

The hrst main result of this paper is on the global existence, uniqueness and well- 
posedness of strong solutions to the Cauchy problem of system (ll.5l) - (ll.9l) : 

Theorem 1.1. Suppose that M.lCh) holds, and the initial data 


(Mo,no,^e,o,ge,o) e with V ■ Mq = 0. 


( 1 . 11 ) 


Then, we have the following: 

(i) There is a unique global strong solution {u,v,Te,qe) to system U.5\) - I[T7^) . with 
initial data {uo,Vo,Tefi,qefl), such that 
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for any positive time T, here and what follows, we use C{- ■ ■) to denote a general 
positive constant depending only on the guantities in the parenthesis. 

(a) Suppose, in addition to that ge,o < 0, a.e. on then 

sup f Wdtq^Wldt < C {a,Q,T,\\{uo,Vo,T^^o,qefl)\\m) , 

0<t<T £ Jo 


for any positive time T■ 

(Hi) Suppose, in addition to that (VTe^o, Vge,o) ^ L”^(R^), for some m G 

(2, oo), then the following estimate holds 


sup 

o<t<r 


\\{VTe,Vqe){t)\\l, < C {a,Q,T, \\{Uo,Vo,Tep,qe,o)\\Hh ||(VTe,o, Vge,o)||m) , 


for any positive time T, and the unique strong solution {u,v,Te,qe) depends contin¬ 
uously on the initial data, on any finite interval of time. 


Formally, by taking the relaxation limit, as e —)■ O’*", system fll.5p - fll.9|) will converge 


to the following limiting system 

dtu + {u ■ V)m — fiAu + Vp + V • (n ® n) = 0, (1-12) 

V-M = 0, (1.13) 

dtv + (u ■ V)v — piAv + (v ■ V)u = — ^ —V(Te — qe), (1-14) 

1 + a 

dtTe-\-u-VTe - {I - Q)'V ■ V = 0, (1-15) 

dtqe + u ■ Vge + (Q + «)V • n < 0, (1-16) 

qe < 0, (1.17) 

dtqe + u ■ Vge + {Q + «)V • n = 0, a.e. on {qe < 0}. (1-18) 


Note that equation fll.9p is now replaced by three inequalities fll.l6l) - fll.l8|) . 

Inequality fll.lbp comes from equation fll.91) . by noticing the negativity of the term 
— h^g+, while inequality fll.l7p is derived by multiplying both sides of equation 
fll.9p by e, and taking the formal limit £ —)■ O’*". Inequality fll.lSp can be derived by 
the following heuristic argument: Let {us,Ve,Tee,qee) be a solution to system fll.Sp - 
fll.9p . and suppose that {ue,Ve,Tee,qee) converges to {u,v,Te,qe), with qe < 0; for 
any compact subset K of the set {{x,y,t) G R^ x (0, oo) | qe{x,y,t) < 0}, since qee 
converges to qe, one may have qee < 0 on 77, for sufficiently small positive e; therefore, 
by equation fll.9l) . it follows that dtqee + ■ V^ee + (Q + «) V • = 0, a.e. on K, from 

which, by taking £ ^ O"*", one can see that (ll.lSp is satished, a.e. on K, and further 
a.e. on {qe < 0}. 

The other aim of this paper is to prove the global existence and uniqueness of 
strong solutions to the limiting system fll.12p - fll.18p . and rigorously justify the above 
formal convergences, as e ^ O’*". Strong solutions to system fll. 1211 - 01.181) are dehned 
in the similar way as those to system fll.5p - fll.9l) . 
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Theorem 1.2. Suppose that holds, and the initial data 

{uo,vo,Tefi,qep) e V-Mo = 0, qe,o < 0, a.e. onR^. (1.19) 

Then, there is a unique global strong solution {u,v,Te,qe) to system 
with initial data (mq, fo, ^e,o, 5'e,o); such that 

sup \\{u,v,Te,qe)it)\\Hi + [ {\\{u, v)\\h 2 + \\Vu\\oo + \\{dtU, dtV, OtT^, dtqe)\\l) dt 
0<t<T Jo 

<C {a,Q,T, |l(Mo,^^0,^e,0,ge.o)||Ri) , 

for any positive time T■ 

If we assume, in addition, that (VTe^o, Vge.o) £ for some m G (2, oo), 

then we have further that 

sup ||(VTe, Vge)(t)||m < C {a,Q,T, |1(Mo,-^ 0 , ^e, 0 , ge,o)||Rl, ||(VTe,o, Vge,o)||m) , 

0<t<T 

for any positive time T, and the unique strong solution {u,v,Te,qe) depends contin¬ 
uously on the initial data. 


Theorem 1.3. Suppose that M.lfh) holds and the initial data 
(uo, Vq, Tefl, qe,o) G H^{R^), V ■ Mo = 0, 

(VTe,o, Vge,o) e qefl < 0, a.e. on 

for some m G (2, oo). Denote by {us,Vs,Tes,qee) one? {u,v,Te,qf) the unique global 
strong solutions to systems / (j.5]) -/ l77PI) and respectively, with the same 

initial data {uq, Vq, q^fl ). 

Then, we have the estimate 

sup II (Ug - U,Ve-V,T^e - Te,qee “ (le)it)\\l 

0<t<T 

+ £ (||(V(m, - m), V(m, - m))||2 + dt < Ce, 

for any finite positive time T, where C is a positive constant depending only on 
a,Q,m,T, and the initial norm || (mq, "Oo, ?e,o, 7;,o)||ni + ||(VTe,o, Vge,o)||m- 
Therefore, in particular, we have the convergences 

{Ue, V,) (n, v) m L°°(0, T; L‘^{R‘^)) n L‘^{0, T; H^{R‘^)), 

{Tee, Qee) ^ {Te, ^e) ^n L-(0, T, L\R'^)) , qfe ^ 0 m L\0, T, L\R‘^)) , 
for any positive time T, and the convergence rate is of order 0{y/e). 


Remark 1.1. (i) In the absence of the barotropic mode, global existence and unique¬ 
ness of strong solutions to the inviscid limiting system was proved in ISw . and the 
relaxation limit, as e O’*", was also studied there, but the convergence rate was 
not achieved. Note that in the absence of the barotropic mode, the limiting system is 
linear, while in the presence of the barotropic mode, the limiting system is nonlinear. 
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(a) Existence and uniqueness of solutions to the limiting system 
without viscosity, was proposed as an open problem in E3/, and also in \3a ] . 

Notably, Theorem \l.^ settles this open problem for the VISCOUS version of 
U.18\) . Note that we only add viscosity to the velocity equations, and we do not use 
any diffusivity in the temperature and moisture equations. 


Remark 1.2. Global well-posedness of strong solutions to a coupled system of the 
primitive equations with moisture (therefore, it is a different system from those con¬ 
sidered in this paper) was recently addressed in where the system under consider¬ 
ation has full dissipation in all dynamical equations, and in particular has diffusivity 
in the temperature and moisture equations. Note that we do not need any diffusivity 
in the temperature and moisture equations in order to establish global regularity of 
the systems considered in this paper. It is worth mentioning that the global regularity 
of the coupled three-dimensional primitive equations with moisture and with partial 
dissipation is a subject of a forthcoming paper. 

The rest of this paper is organized as follows: in section |2l we state and prove several 
preliminary lemmas, while the proofs of Theorem 11.11 Theorem 11.21 and Theorem 11.31 
are given in section |3l section H] and section |5l respectively. The last section is an 
appendix in which we prove some parabolic estimates that are used in this paper, 
and which are of general interest on their own. 


2. Preliminaries 


We will frequently use the following Ladyzhenskaya inequality (see, e.g., [26 


ll/IIW)<C^II/ll 1,2 (R2) liv/ll L2(R2)) V/ G H\R^). 

The following lemma on the Gronwall type inequality will be used to establish the 
global in time a priori estimates to the strong solutions to system fll.5p - fll.9l) later. 


Lemma 2.1. Given a positive time T, a positive integer n and positive numbers 
ri G [1, oo), 1 < i < n. Let Oq, ai and h, 1 < i < n, be nonnegative functions, such 
that Oo, tti G L°°((0, T)) and bi G T^((0, T)). Suppose that the nonnegative measurable 
function f satisfies 

f{t) < ao{t) + Vai(t) f f bi{s)ffis)ds 
i=l V./0 

for any t E [0, 7~]. Then, the following holds 



||/IU-((o,r)) <(i7 + l)'’ iaolLexp 


(r, + l)’-5^||a.||;,(l + ||6.||ir+' 

i=l 




where r = maxi<i<„ri, and || ■ ||i and || • ||oo denote the L^((0,T)) and L°°((0,T)) 
norms, respectively. 
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Proof. By the Holder and Young inequalities, we deduce 


biis)p{.s)ds 




Uo 


{s) [bl{s)f{s)] ds 


< 


bi{s)ds 


bi{s)f{s)ds 


< (1 + ll^illi) bi{s)f{s)ds] , 
for 1 < i < n. Therefore, by assumption, we have 

n 

/(t) <||ao||oo + 5]] ||ai|loo(l + ll&illi) ( / bi{s)f{s)ds] , 


i=l 


from which, taking the r-th powers to both sides of the above inequality, and using the 
elementary inequality (X]r=o — (^ + 1)^”^ ^ 27=0 where q are positive numbers, 
we arrive at 


rm <(n+ l)’'-‘||aollL + (" + 1)’'"' E 

i=l 


bi{s)f{s)ds 


Applying the Gronwall inequality to the above inequality, we have 


r{t)pn + iy iaolLexp (n +1)^ ^ l|ai|L(l + / bi{s)ds 

I i=i 

<(n + l)'-iao||;;oexp |(n + 1)^"^ ^ l|ai|L(l + 


\r+l 


2 = 1 


from which, taking the r-th power root to both sides of the above inequality, and 
taking the supremum with respective to t over (0, T), one obtains the conclusion. □ 


The next lemma will be employed to prove the uniqueness of strong solutions. 


Lemma 2.2. Given a positive time T, and let mi, m 2 and S be nonnegative functions 
on (0, T), such that 

mi. S'G L^((0, T)), m 2 G L^((0, T)), and S'> 0, a.e.on{0,T). 


Suppose that f and G are two nonnegative functions on (0, T), with f being absolutely 
continuous on [0,T), and satisfy 


1 fit) + Gif) < mi{t)f{t) + m 2 {f) 

f{t)G{t) log+ (^|||) 

\ /(O) = 0, 



a.e. on (0, T), 
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where log’*'= max{0, log^;}, for z G (0, cxd), and when G(t) = 0, at some time 
t G [0,T), we adopt the following natural convention 

(H)= }% " (^)=“■ 

Then, we have f = 0 on [0, T). 


Proof. Suppose, by contradiction, that there is some time G (0,T), such that 
/(t*) > 0. Recalling that / is absolutely continuous on [0,T), by the property of 
continuous functions, there must be a time to G [0,f*), such that /(to) = 0 and 
fit) > 0, for any t G (to,t*]. In the rest of the proof, we will focus on the time 
interval [to,t*). For any a G (0, cxd), one can easily check that 

log"*" < —, for 2 : G (0, cxd). 
ae 

Recall the Young inequality of the form ab < ^ for any nonnegative numbers 
a, b, and for any p,q E (1, cxd), with ^ ^ = 1. Thanks to the above inequality, and 

choosing a G (0,1), it follows from the assnmption and the Yonng inequality that 


f’ + G < 

< 


< 


mif + m 2 


fG 


S 

ea \G 


a 1 — <y 

= mi/ + m2S^G^~ 




f \2 

m2^^ ( — 


ea 


1 + 0 - 


1-a^ , 1 + a / / ^ '+" 

2 


mi/ + mo+"5'i+" ( — 


ea 


mif + G + m 2 ^‘"3^+^^ ( , a.e. on(0,T). 


a 



1 

2 


Note that the arguments used in the above ineqnality are for the time when G{t) > 0; 
however, for the time when G{t) = 0, recalling that we nnderstood the term involving 
G as zero, therefore, the above ineqnality result holds trivially. Therefore, we obtain 

/'<mi/ + mr^^^ ^ , 


for any a G (0,1), and for a.e. t G [to;fo)- Recall that fit) > 0, for t G (to,R)- 
Dividing both sides of the above ineqnality by /i+^, then one can dednce 



ai+‘^ t/- ^-2_ 

< ——mi/i+- + ——m2+"S^i+'^ 

1 + a 1 + a 

(T 'T (T ^ cr 

< -mif“ + a“mn+'^S'“, 

“ 1 + a ^ 
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for a.e. t G (to,t*). Applying the Gronwall inequality to the above inequality, and 
recalling that /(to) = 0, it follows from the Holder inequality that 


/i+'^(t) < ae 1+'^-^*0/ 7712 ^'" {s)S^+<^ {s)ds 


'to 


< ae 1+'^-^*0 I / ml{s)ds 

'to 


1 

1 + <T 


S{s)ds 


1 + <T 


from which, taking the ^^-th power to both sides of the above inequality, one obtains 

/(t) < (Te-^*o ml{s)ds^ j S{s)ds, 

for any t G [to,t*), and for any a G (0,1). Recall that m 2 G L^((0,T)), by the 
absolute continuity of the integrals, there is a positive number 17 < t* — to, such that 
m\{s)ds < 1, for any t G [to, to + rj). Therefore, the above inequality implies 


/(t) / S{s)ds, 


for any t G [to, to + 17 ), and for any a G (0, 1). By taking a —)■ O'*', this implies that 
/ = 0, for any t G [to, to + 17 ), which contradicts the assumption that /(t) > 0, for 
any t G (to,t*). This contradiction implies that there is no such t* G (0,T) that 
/(t*) > 0 , in other words, recalling that / is a nonnegative function, we have / = 0 
on [0,T). This completes the proof. □ 


We also will use the following elementary lemma. 

Lemma 2.3. Let LL he a measurable set of positive measure, and f be a mea¬ 

surable funetion defined on H. Suppose that, for any positive number rj, there is a 
measurable subset Erj of Q, with \Ej^\ < rj, sueh that / = 0, a.e. on Q\ Ej^. Then, 
/ = 0, a.e. on H. 

Proof. Suppose, by contradiction, that the conclusion does not hold. Then there is 
a subset E of H, with 0 < |i?| < 00 , such that |/| > 0 on i?, here \E\ denotes the 
LApebessgue measure of the subset E. Then, for 77 = -1^, by assumption, there is a 
subset Er, of H, with \Ef\ < 77 , such that / = 0 on H \ This implies that E C E^, 
and thus 

IT'I < \EP\ <11= 

Therefore, \E\ = 0, which contradicts the assumption that \E\ > 0. This contradic¬ 
tion implies the conclusion of the lemma. □ 
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3. Global existence and uniqueness of the system with positive £ 

In this section, we will prove the global existence and uniqueness of strong solutions 
to the Cauchy problem of system fll.5l) - fll.9|) . for any positive e. Several ^-independent 
a priori estimates will also be obtained. 

Let’s start with the following result on the local existence and uniqueness of strong 
solutions to the Cauchy problem to system fll.5l) - fll.9p . 

Proposition 3.1. Suppose that M.1(A) holds. Then, for any initial data 

{uo, vq, Te,o, qe,o) € with V ■ mq = 0 , 

there is a unique local strong solution {u,v,Tf.,qf) to system / lj.5)) -/ f77P)) . on x 
(0,T), with initial data (mq, uo, 7’e,o, <?e,o); where the existence time T depends on a, 
Q, e and the initial norm ||(mo, uo, 7’e,o, <?e,o)||Hi ■ 

Proof, (i) The existence. The existence of strong solutions to system fll.5p - fll.9l) . 
with initial data (mq, tq, Te,o, qe,o) can be proven by the standard regularization argu¬ 
ment as follows: (i) adding the diffusivity terms —rjATe and —rjAqe to the left-hand 
sides of equations fll.Sp and fll.9p . respectively, in other words, we consider the fol¬ 
lowing regularized system 

dtu + {u ■ V)m — Au -I- Vp -I- V • (u ® n) = 0, 

V • M = 0, 

< dtv + {u ■V)v - Av + {v ■V)u = j^V{Te - qe), (3.1) 

dtTe + u ■ VTe - (1 - Q)V • V - rjATe = 0, 
dtqe + U ■ Vqe + (Q + a)v • v - pAqe = -^qf ; 

(ii) for each p > 0, the Cauchy problem of the regularized system fl3.1l) . with initial 
data {uo,vo,Tefl,qe,o), has a unique short time strong solution {u^'^\v^'^\Te^\q^e'^), 
which satisfies some 77 -independent a priori estimates, on some 77 -independent time 
interval ( 0 , T), for a positive time T depending only on on a, Q, e and the initial norm 
\\{uo,Vo,Tefl,qe,o)\\m] (hi) thanks to these 77 -independent estimates, by adopting the 
Cantor diagonal argument, one can apply the Aubin-Lions lemma and take the limit 
77 —)■ O’*" to show the local existence of strong solutions to the Cauchy problem of 
system fll.5p - fll.9p . with initial data (mq, Vq, Tep, qe,o)- Since the proof is standard, we 
omit it here; however, the key part of the proof, i.e., the relevant a priori estimates, 
are essentially contained in the ’’formal” proofs of Propositions 13.2113.51 below. As 
it was mentioned above, these formal estimates can be rigorously justified by estab¬ 
lishing them first, to be 77 —independent, for the regularized system fl3.ip and then 
passing with the limit as 77 —)• 0 ^. 

(ii) The uniqueness. Let (m, v, Te, qe) and (h, v, Te, (je) be two strong solutions to 
system fll.5p - fll.9p . with the same initial data (mq, tq, <?e,o), on the time interval 
(0,T). Define the new functions 

{6u,5v,6Te,6qe) = {u,V,Te,qe) - {u,V,fe,qe). 
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Then, one can easily check that 

dt5u + {u ■ V)5u + {5u ■ V)u — A6u + Vdp + V ■ {v ^ 6v + 6v ^ v) = 0, (3.2) 

V ■ (5m = 0, (3.3) 

dt6v + {u ■ 'V)6v + {6u ■ V)h — A6v + {v ■ 'V)6u 

+{6vV)u = j^V{6T,-6q,), (3.4) 

dt6T, + u ■ V6Te + 6u ■ VT; - (1 - Q) V • (5m = 0, (3.5) 

dt^Qe + u ■ V5qe + (5m • Vge + (Q + «)V • (5m = -^{qt -qt)- (3-6) 

Since equations fl3.2l) - fl3.5p hold in L^( 0 , T; T^(R^)), we multiply equations fl3.2l) . 
(13.4p and (13.5p by (5m, 5 v and ^Tg, respectively, and integrating over R^, then it follows 
from integration by parts that 

1 d 


2 dt 


[((5m • V)m + V • (m <8) ^m + 5m (g) m)] • Sudxdy 


[(5m ■ V)m + (m ■ V)5m + (5m • V)m] • 5m + ^ dxdy 


1 -|- 


’} 


- / [5m ■ Vfg5Te - (1 - Q)V ■ 5M]5Tgdxdi/=: J. 
By the Young inequality, we deduce 


/< / [|5m||Vm| + (|m| + |m|)|V5m| + (|Vm| + |VM|)|5M|]|5M|(ia;(ii/ 


[|5m|(|Vm| + |Vm|) + |m||V5m|]|5m| + 


|V5m| 
1 T Q( 


(|5Te| + |5ge|) > dxdy 


+ f [\6u\\Vf^\\6T^\ + {l-Q)\V5v\\6Te\]dxdy 
<11 i\V6u\^+ \V6v\^)dxdy + C [ [(|Vm| + |Vm| + |Vm| + |m 

2 Jr2 J^2 

+ |mP)(|5mP + |5mP) + |5Tep + |5geP + \'^Te\\Su\\5Te\]dxdy, 


and thus 

YC [ [(|Vm| + |Vm| + |Vm| + |Mp + |hn(|5M|2 + |5Mn 
+ |5Te|2 + |5gg|2 + \Vfe\\6u\\6Te\]dxdy. 


(3.7) 
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Multiplying equation (13.6^ by 6qe, integrating the resultant over then it follows 
from integration by parts and the Young inequality that 

[ {(it - Qt)iQe - qe)dxdy 
z at e j ^2 

= — {6u ■ Vqe + (a + Q)V ■ 6v)6qedxdy 

<]\\^dv\\l + C f {\6qe\^ + \Vqe\\Su\\6qe\)dxdy, 

4 JR2 

from which, noticing that the function is nondecreasing in thus — q'^){qe — 
qe) > 0, and one obtains 

^WdQeWl < ^||V5u||2 + C f (|5geP + \Vqe\\6u\\6qe\)dxdy. 

at z j ^2 

Summing the above inequality with fl3.7|) yields 

< c[ [(|Vh| + |Vh| + |Vu| + |u |2 + |un(|5u|2 + |5t;|2) 

+ \5T,\^ + \5q,\^ + \Vfe\\Su\\6T,\ + \Vqe\\Su\\5qe\]dxdy, (3.8) 

from which, by the Holder, Ladyzhenskay and Young inequalities, we deduce 

||| {6u, 6v, 5Te, 6qe)\\l + ^(11 V5u||^ + || V5u||2) 

< C(||(Vh,Vh,Vu)||2 + ||(u,h)||2)||(5u,5u)||2 

+C'||(5Te,5ge)||^ + C'||(Vfe,Vge)||2||5n||oo||(^Te,^ge)||2 

< C(||(Vh,Vh,Vu)||2 + ||(u,u)||2)||(5u,Hll2||(V5u,VHIl2 

+C||(5Te,5ge)||^ + C^||(Vfe,Vge)||2||5n||oo||(^Te,5ge)||2 

< ^||(V5u,VMIl2 + C^(ll(Vh,Vu,Vu)||2 + ||(h,h)||^)||(5u,5u)ll2 

+q|(5Te,5ge)||^ + qi(Vfe,Vgg||2||(5n,Hlloo||(5Te,5g,)||2. 

Therefore, one has 

^\\{Su,Sv,ST„Sqc)\\l + i||(i!i,iw)|ll,i 

< C(1 + ||(fi,«)||l+ ||(Vfi.VS,Vt.)ll8 \\{Su,6v,ST„SqMl 

+c||(vL,vaih||(fe.fo)IUI(i 5 r.,i 9 ,)|| 2 . 


(3.9) 
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Recalling the following Brezis-Gallouet-Wainger inequality (see 00) 


||/||l°°(k 2 ) < C'||/||j^i(r 2 ) log 2 r|K !|—-—- + e 


1 / ll/l|g2(R2) \ 

Vll/llw 


and denoting U = {u,v),U = {u,v) and 6U = {6u,6v), we have 




(3.10) 


where 


S{t) — \\U\\h'^ + \\U\\h^ + e(||t/||jyi + \\U\\Hi)- 


Note that, when 6U = 0, fl3.10p still holds, as long as we understand the quantity on 
the right-hand side as zero, in the natural way as in Lemma 12.21 


Denoting 



mi = G (1 + ||(h,h)||^ + ||(Vh, Vh, Vn)||2) , = C\\{Vt,Vm2, 

then it follows from fl3.9l) and fl3.10p that 



Here, at the time when G{t) = 0, the term involving G{t) on the right-hand side 
of the above inequality is understood as zero, as it was in Lemma 12.21 Recalling 
the regularities of {u,v,Te,qe) and {u,v,Te,qe), one can easily check, thanks to the 
Ladyzhanskaya inequality, that mi,S G L^((0,T)) and m 2 G L^((0,T)). Therefore, 
we can apply Lemma [2.21 to conclude that / = 0, which proves the uniqueness. □ 

For the rest of this section, we always suppose that (m, v , T^, qe) is the unique 
strong solution to system fll.5p - fll.9p . on x (0, T), for some positive time T, with 
initial data (mq, vq, Tg^o, qe,o)- We are going to establish several ^-independent a priori 
estimates on {u,v,Tg,qg). Before performing these a priori estimates, we point out, 
again, that the arguments being used in the proofs of Propositions I3.3H3.51 below, 
are somewhat formal, because (u, u, Te, qg) may not have the required smoothness for 
justifying the arguments. However, one can follow the same arguments presented 
in the proofs of Propositions I3.3H3.5I to establish the same a priori estimates to the 
regularized system fl3.ip . for which the solutions fulhll the required smoothness, and 
then take the limit —)■ O’*", recalling the weakly lower semi-continuity of the relevant 
norms, to obtain the desired a priori estimates on (u, u, Te, qg). 

Let’s start with the basic energy equality stated in the following proposition. We 
observe that here we have energy equality, instead of inequality, as in the case of 
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strong solutions of the Navier-Stokes equations. Observe, however, that for the rest 
of the proof of the main result it is sufficient to have energy inequality. 


Proposition 3.2. We have the following estimate 

|U||2 I ||y||2 I W'^eWl _ WqeWl \ 

2dty "2^ (i + a)(l-Q) ^ (l + a)(Q + a)y 

II II2 

+ \\Vu\\l + \\Vv\\l + ^^2^ = 0, 

e{Q + a) 

for any t G (0, T). 


Proof. Multiplying equations fll.Sp and fll.7p by u and t, respectively, summing the 
resultants up and integrating over then it follows from integration by parts that 

l^i\\u\\l + \\v\\l)+ \\Vu\\l + \\Vv\\l =f {qe-T^)V-vdxdy, (3.11) 
2dt 1 + a Jr2 

where we have used the following fact 


/ [V • (n 0 n) • M + (t • V)m • v]dxdy 

JR2 



(t 0 t) : Vu\dxdy = 0. 


Multiplying equation fll.Sp by (1 + a) ^(1 — Q) and integrating over M^, then it 
follows from integration by parts that 


d 


2 (l + a)(l-g) dt 


IITII^ - 


1 “f- 


TeV * vdxdy = 0. 


(3.12) 


Multiplying equation (II. 9p by (1 + a) ^{Q + a) ^Qe, and integrating over M^, then it 
follows from integration by parts that 


d 


2(1 + cr)(g + o) dt 


h 


2 

e|l2' 


1 "h O 


-vdxdy = 


siQ + o) 


\qf\^dxdy. (3.13) 


Summing fl3.11l) - fl3.13|) up yields the conclusion. 


□ 


As an intermediate step to obtain the L°°(0, T; i7^(M^)) estimate for {u,v,Te,qe), 
we prove the L°°{0, T; estimate in the next proposition. 

Proposition 3.3. Denote U = {u,v). Then, we have the estimate 

sup \mT„q,)\\t+ r(||T|VC/||0||V.i,||l)*<C, 

o<t<r Jo ^ 2 

for a positive constant C depending only on the parameters a, Q, T and the initial 
norm || (f^o, Te^oi ?e,o)||L2(R2)nL4(R2), and in particular, C is independent of e. 
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Proof. Multiplying equations (ll.5p and (II.7p by \U\‘^u and \U\‘^v, respectively, sum¬ 
ming the resultants up and integrating over then it follows from integration by 
parts and the Holder inequality that 


1 d 
4 dt 


1 + 


\U\^\VU\^ + -\V{\Un^]dxdy 


= / [pV ■ {\U\\) - (V ■ (u 0 u)) • \U\\ 

rri 

— {v ■ W)u ■ \U\^v -f ^■ {\U\^v)]dxdy 
1 + a 

<3 / [H|HnVH| + |HnVH| + (|g,| + |Te|)|HnVH|]da:di/ 

./R2 

<3(||p|| 4+ |||f/nL+ ||Te||4+ |ke||4)||t/||4|||f/|VH||^. 


(3.14) 


Applying the divergence operator to equation fll.Sp . in view of fll.bp . one can see that 


—Ap = V -V ■ {u ® uV ® v). 


Note that p is uniquely determined by the above elliptic equation by assuming that 
p —)■ 0, as (x, y) —)■ oo. Thus, by the elliptic estimates, one has 

IIpIU < C\\u ^ u + V ^ v\\4 < (7111171^11^. 

Substituting this estimate into fl3.14p . and using the Ladyzhenskaya and Young in¬ 
equalities, one deduces 

~l|C'll3+||IDvy|0i||v|c/ni^ 

<c(|||£/p|L + lir.iu + ||94|4)||t;||4|||f/|V!/||, 

<C-(|||(7|2|||||V|(7|2||j + IIT.IU + |l94l4)l|f/||4|||f/|Vf/||, 

+ cillf/llJdir.ii^ + \\q,f,) + ||(7|ia 
<C-(1 + ||C/||1)(||U||^ +II,.11^+1117111) 

+|(iTi^c/|iYi7if^ni7 

and thus 

7||f/||l + 2\\\U\VU\\l < C(1 + ||(7||1)(||T.||1 + II,.111 + ||f/||l). (3.15) 

Multiplying equation fll.Sp by \Tff‘Te, and integrating over R^, then it follows from 
integration by parts and the Holder inequality that 

1 7 n 

--||T,||^ = (1-Q) / V-u|Te|2T,dxdp<(l-Q)||Vu||4||Te||i 
4 at 7]r2 
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which implies 


^||T,||2<2(1-Q)||Vt||4||T,||4. 


Similar manipulation to equation fll.9p yields 


A 

dt 


Ik.llS < 2((3 + a)||Vt>||4||?.||4. 


(3.16) 

(3.17) 


Summing fl3.15l) - fl3.17l) up, and integrating the resultant in t yield 

(ITIIi+IIT||l+|l94||5)(i) + 2 f \\\U\VU\\lds 

Jo 

— ll^o||4 + ll^e,o||4 + lke,o||4 + 2(1 + «) / ||Vu|| 4 (|| Tg || 4 + ||ge|| 4 )ds 


+ / il+\\U\\t)m\\t+\\Te\\l+\\qe\\l)d^ 


(3.18) 


for t E [0, T). 

We need to estimate the term f* || Vn|| 4 (||Te ||4 + ||( 3 'e|| 4 )ds on the right-hand side of 
(13.181) . To this end, applying Lemma [62] (in the Appendix section) to equation (11.71) 
yields 

2 


\\Vv\\lds <C 


l|Vno|l2 + 


|17|V17||> + / {\m + \\qe\\t)ds 


, (3.19) 


for all t G [0,T), where C is a positive constant independent of t. Thanks to this 
estimate, it follows from the Holder and Young inequalities that 

2{1 + a) [ ||Vn||4(||Te||4-|-IlgelU)^^ 




llVnll^ds 


m\\l+\\q4l)ds 


<Ct -2 ||Vno||^+ / \\\U\VU\\Us 


{\\Ut + \\qe\\t)ds 


+ Ct-^ I / (||Tg||^+||gg||^)ds 


< 


iHiVt/ll'ds + C'l / i\\Te\\t + \\qe\\t)ds] +C'||Vno|| 


1^ 


Substituting (I3.20p into (I3.18p . and denoting 

/(() = (ITII* + liTjiJ + ||«4I5)(() + r |||(7|vc/||"<ii,, 


(3.20) 
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we have 

/(i) < /(O) + C||+ C(j‘ f(s)ds'\‘ +cj\l + ||t;||S/(s)<is. 

for all t G [0,T). By Proposition 13.21 and using the Ladyzhenskaya inequality, one 
can easily check that /q^( 1 + ||17||4)(it < C, for a positive constant C depending only 
on a,Q,T and the initial norm || (mq, uq, Te,o, 5 'e,o)|| 2 - Therefore, applying Lemma [2T] 
to the above inequality, one obtains 

pT 

sup \\{U,T,,q,){t)\\l+ / \\\U\VU\\\dt<C, 

0 <t<T Jo 

and further, recalling (Id.lbh . proves the conclusion. □ 

Thanks to the a priori estimate stated in the above proposition, one can immedi¬ 
ately obtain the L°°(0, T; estimate on u as stated in the following proposi¬ 

tion. 

Proposition 3.4. ITe have the following estimates 

sup ||VM(t)|| 2 + [ \\Au\\lds<C, 

0 <t<T Jo 

for a positive constant C depending only on the parameters a, Q, T and the initial 
norm ||Mo||iri(]R 2 ) + || ('I'o, Te^o, ?e,o)||L2(]R2)ni4(R2), and in particular is independent of e. 

Proof. Multiplying equation (II.5p by —Am, and integrating over M^, then it follows 
from integration by parts that 

a 4 ||Vm ||2 + ||Am ||2 = [ [{u-V)u + V{v(^v)]-Audxdy 
i at j^2 

<3 [ |t/||V17||AM|da:di/<i||AM||2 + C|||f/|V17|||2, 
where, again, U = {u,v), and thus 

L||v«||:i + ||A«||:i<c|||(7|vf/|||j, 

for all t G [0, T). From which, in view of Proposition 13.31 the conclusion follows. □ 

Finally, we are ready to prove the L°°(0, T; estimate on (n, Tg, qe), that is 

the following proposition. 

Proposition 3.5. The following estimate holds 

sup ||(VM,VTe,Vge)(t)||^+ r (\\Av\\l+^^^^ + \\Vuu)dt<C, 

0<t<T Jo \ ^ J 
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where C is a positive constant depending only on a, Q, T and the initial norms 
\\{uo,Vo,Tefl,qep)\\m, and in particular is independent of e. 


Proof. Multiplying equation fll.Sj) by — ATg, and integrating over then it follows 
from integration by parts and the Holder inequality that 

\-^\\S7Te\\l ={l - Q) [ WTe-Wiy ■ v)dxdy - f diU ■ WTediTedxdy 

^ dt 7k2 7r2 

<(1 - Q)||An||2||VTe||2 + ||Vn||oo||VT,||2. 

Similarly, one can derive from equation fll.Qp that 

^^11 WII 2 + ^11 vy 11 ^ < (a + |5)I|AulhllVfcih + ||V«|U||V,.||^. 

Summing the previous two inequalities up yields 

~(l|VT.||^ + llv«.||=) + i^llvyii^ 

<(1 + a)|| Awlhdl VT.Ib + llVg.lh) + ||Vti|U(|| VT.II" + II ViJ.I© 

<i|| A't-ili + lliv«iu + 2(a + i)T|| vr.iii + II V«.|ll), 

and thus 

^(livr.||^ + ||v,.||i) + i^l|v,.+||^ 

< 2|||V«|U + 2(a+l)Tl|VTd|I+||V*|l^) + i||AK||^. (3.21) 

Multiplying equation fll.7p by —An, and integrating over R^, then it follows from 
integration by parts and the Young inequality that 


l±\\Vv\\l+\\Avg = 


1 -|- 


V(Te - Qe) - {u ■ V)n - (n ■ V)m 


Avdxdy 


1 , 


<-||A!.||^ + c IIVTdIl + IIV«dl^ + IIUIVC/ 


and thus 


fL 

dt 


IIVHI 2 + ||At’ll^ < c (||vrdl 2 + l|V<;.||^ + ||U|VC/||^), (3.22) 


Summing fl3.2ip with fl3.22p up yields 


|||(Vn, VT„ Vge)||^ + ||An||2 + ^\\Vqt\\l 


< C'll|t/|Vt /||2 + CdlVnlloo + 1)(||VTell^ + IIVg, 


ell^), 
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from which, by the Gronwall inequality, and using Proposition 13. 3 [ one obtains 

'T' / 

sup ||(Vi,,Vr.,VgOWIl 2 + / (||A<’II 2 + —lIVyilSMi 


0<t<T 

rT, 


e 

rT 


(||v.|U+i)dt ||(Vt;o,VTe,o,Vge,o )||2 + C' / \\\U\VU\\ldt 


<C'(a,Q,T, |l(t/o,re,o,ge.o)||Hi)exp < G / (|| Vm||oo + l)dt ^. 


rT 


(3.23) 


To complete the proof, one still need to estimate ||Vm||oo<^^- If follows from 
Propositions I3.3H3.4I and the Ladyzhenskaya inequality that 


f-r 


r-r 


/ {\\Vu\\t+\\Vv\\t)dt<C iWVuUWAuU + \\Vv\\t)dt 
Jo Jo 

<C{a, Q, T , ||(mo, Tefl-, qe,o)\\m)- (3.24) 

We decompose u as u = it + it, where it and u, respectively, are the unique solutions 
to the following two systems 

dtu — Ait + Vp = —{u ■ V)m — V ■ (n ® n), 

V-'u = 0, (3.25) 

u\t=o = 0, 

and 

dtii — Ait + Vp = 0, 

V-M = 0, (3.26) 

it\t=o = uq. 

We are going to estimate it and it. Let’s first estimate it. By the L'^(0, T; W^’^) 
type estimates for the Stokes equations (see, e.g., Solonnikov |39|, l40|), we have 

IKcIih, Ah)||L5(]R2x(o,r)) A G|||f/|Vt/||i9(]R2x(o,r))) 
for any q G (l,oo), and thus it follows from the Holder inequality and Gagliardo- 
Nirenberg inequality, ||(p ||?2 A Cllv^lllllV(p|| 4 , fl3.24p and Proposition 13.31 that 


rT 


fT 


-r 


||Afi||^A<C/ |||C/|V£;||p«<C / llVt/IIJIIt/ll^jdi 

Jo Jo 

rT 


<C 


II VGlI^llGlI^dt < G(a, Q, T, |1(mo, Te,o, qe,o)\\m). 


One can deduce easily from equation fl3.25|) . by using Proposition 13.31 that 


-r 


rT 


sup ||VM(t)|| 2 + / \\Au\\ldt<C / |||G|Vf/||^dt 

0 <t<T Jo Jo 

<C{a, Q, T, II ( mo , Vo, Te,o, qe,o)\\H^). 
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Thanks to the above two estimates, it follows from the Gagliardo-Nirenberg, 

1 3 

(711 v?II 2 IIA(^III, and the Holder ineqnalities that 


< 


r-T 


-r 


\\Vu\\oodt<C \\VuM\\Au\\idt 


<C 


fT 


WVuWldt 


rT 


II AnIIgdt) T5 


<C{a, Q, T, II (mo, uo, Te,o, ge,o)||Hi)- 


(3.27) 


Next, we estimate u. Mnltiplying eqnation fl3.26p by (tA^ — A)u, and integrating 
the resnltant over R^, then it follows from integration by parts that 

^^(l|Vh||| + \\ViAu\\l) + ^\\Au\\l + \\ViVAu\\l = 0. 


Therefore, we have 


rT 


snp 

0<t<T 


(livens+ ||v^Afl||^)+ / (||Afi||^ + ||AVAfi||l)*<||V«„||2. 


Thanks to this estimate, it follows from the Gagliardo-Nirenberg (Agmon), 
1 1 

G||(p||| II A(p|||, and Holder ineqnalities that 


< 


r-T 


\\Vu\\^dt<C [ ||Vh|||||VAM|||dt = G [ ||VM|||||\/tVAh|||r3dt 


<G 


rT 


WVuWldt 


fT 


\\\/tW Au\\\dt 


rT \2 

t ^dt 


'0 


AGr^llVnolllllVnolllr^ = Cr^WVuoh- 

Gombining the above estimate with fl3.27p . one has 


fT 


fT 


||Vn|loodt < / (||Vm||oo + ||Vh||oo)dt < C{a,Q,T, ||(Mo,uo,Te,o,ge,o)||RU 


which, when snbstitnted into 03.231) . yields the conclusion. 


□ 


As a corollary of Propositions 13.2f[331 we have the a priori estimate to (u, n, Te, ge), 
as stated in the following: 

Corollary 3.1. Suppose that U.1(J\) holds, and the initial data 

{uo,vo,Tefi,qefi) E V ■ rto = 0. (3.28) 

Let (u.v.T^.aA be the uniaue strona solution to sustem /TOD-zTOI). on R^x (o,r), 
0 < T < oo, with initial data (mq) ^e, 0 ) <?e,o)- Then, the following hold: 
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(i) We have the estimate 

snp \\{u,v,Te,qe)it)\\jji+ [ + \\{u,v)\\h 2 + ||VM||oo)c^t 

0 <t<T Jo \ ^ J 

+ [ \\{dtu,dtv,dtTe)\\ldt < C {a,Q,T,\\{uo,vo,Tefi,qe,o)\\m) ■ 

Jo 

(ii) Suppose in addition to H3.28^) that q'^^ = 0, a.e. on R^, then we have 

||( 7+||2 pT _ 

sup J—/ \\dtqe\\ldt <C {a,Q,T,\\{uo,vo,Tefi,qe,o)\\m) ■ 
0 <t<T £ Jo 


(Hi) Assume in addition to H3.28\) that (VTe^o, Vge.o) ^ L”*(R^), for some m G 
(2, cxo), then we have the estimate 


sup 

0<t<T 


IKVTe, Vge)||m < C {a,Q,T,m, UK, Uo, Tg^o, ge,o)||j/l, ||(VTe,o, Vge,o)||m) • 


Proof, (i) The estimate on all the terms, except those involving the time deriva¬ 
tives, follow directly from Propositions I3.2H3.51 The desired estimate for {dtu, dtv) 
follows directly from the a priori estimate in Propositions 13.31 and 13.51 by using the 
L‘^{0,T; H^) type estimates to the Stokes and heat equations. By Propositions 13.21 
13.41 and 13.51 it follows from equation fll.Sp and the Sobolev embedding inequalities 
that 

' 7 ~ 7 ~ 

/ \\d,TA\ldt<[ i(i-Q)||VK||l + ||,.||^||vr.|ii|* 

Jo Jo 

T T 

<C + C [ \\u\\l^dt<C + C [ \\ufH 2 dt<C. 

Jo Jo 

(ii) Multiplying equation fll.Qp by dtqe, and integrating over R^, then it follows 
from the Young and Sobolev embedding inequalities and Proposition 13.51 that 

1 (y. d f 

-Tr—^Mt\\l + \\dtqe\\l = - [u-Vqe + {Q + a)V ■v]dtqedxdy 

le at j ^2 

< )l|9Alll + C(||«||Ll|Vfe||l + ||VHIi) 

^ + r '(|| u ||^3 + 1 ), 

from which, by (i), the conclusion in (ii) follows. 

(hi) Applying the operator V to equation fll.81) . multiplying the resultant by 
IVTgp^^VTe, and integrating over R^, then it follows from integration by parts 
and the Holder inequality that 

-^AVT4f)={l-Q) [ \VT,r-^VT,-V{V-v)dxdy 
m dt 
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djU ■ VTplVTf 


I m—2 


diTf,dxdy 


Thus 


<(1 - g)||V2u|U||VTeC"' + ||Vu|UI|VT, 


^lIVTelU < (1 - Q)\\VMm + ||Vu|U||VT, 
dt 


m 
e llm' 


Similarly, one can derive from equation fll.Qp that 


d 

dt 


llV^elU < (a + g)||V^n||m + llVnlloollVge 


Summing the above two inequalities, one obtains 

— (lIVTellm + ||Vge||m) < (1 + <^)||V^n||m + ||Vu||oo(||VTg||m + llVlJellm); 
from which, integrating with respect to t, we have 


||(VTe, Vge)llm(t) < C* / \\V^v\\mds + C || Vu || oo || (VTg, Vge) || (3.29) 

Jo Jo 

for all t G [0, T). 

Applying Lemma EISl see the Appendix section below, to equation fll.7p . and using 
the Sobolev embedding inequality, one deduces 



/ft \ 2 

/ 1 V^n| mds < C 

l|Vt.„|| 2 + / ||(vr„v,.,|u||vi.|,|t.||vi,|)||> 

Jo 

\Jo / 


<(7 


<(7 


l|v-<,l|2+ / (ll(vr.,v,.)||^ + ||'t.|iLl|Vt>|| 2m "h I [i’ll 2m liVnII 

2m) 


II Vn, 


0||2 


(II (VTe, Vge) IIL + II {U, v)\\l, II (Vn, Vn) \\l,)ds 


for any t G [0, T), where C* is a positive constant depending only on m and T, and 
is in particular independent of t G [0,T). By (i), the above inequality implies 

WVMmds <C + c(^j^ IKVTe, VqeWmds 

for any t G [0, T), and for a positive constant C independent of t G [0,T). Substi¬ 
tuting the above estimate into fl3.29l) . and setting f{t) = ||(VTe, Vge)llm(t) yield 


||Vno||2+ / fisYds] + / ||VM||oo/(s)ds 


fit) < c 


for any t G [0, T), where C is a positive constant independent of t G [0, T). Recalling 
(i), and applying Lemma [2.11 the conclusion stated in (iii) follows. □ 
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Now, we are ready to prove the global existence, uniqueness and well-posedness of 
strong solutions to the Cauchy problem of system fll.5|) - fll.9l) : 


Proof of Theorem \l. 1\ The uniqueness of strong solutions follows from Proposi¬ 
tion EU] directly, while the a pnori estimates in (i)~(iii) follow from (i)-(iii) of Corol¬ 
lary [3T1 respectively. Therefore, we still need to prove the global existence of strong 
solutions as stated in (i), and the continuous dependence of the strong solutions on 
the initial date as stated in (iii). 

To prove the global existence of strong solutions, it suffices to extend the local 
solution established in Proposition 13.II to be a global one. By repeating Proposition 
13.11 one can extend the local solution (m, n, Tg, qe) to the maximal interval of existence 
[0, 7Z). Then, we need to show that % = oo. Suppose, by contradiction, that % < oo, 
then we must have 

lim \\{u,v,Te,qe)\\Hi = oo. 


However, by Corollary 13.11 which holds since 7^ < oo, the quantity ||(M,n,Te, qe)\\l. 
is bounded on [0, 7^), which is a contradiction, and thus % = oo. 

We now prove the continuous dependence of the unique strong solutions on the 
initial data as stated in (iii) on any hnite interval [0,7~]. Therefore, we choose ar¬ 
bitrary T G (0,oo), and focus on the interval [0,7~]. Let {u^^\v^^\Te^\qi^'^) and 
be the unique solutions to system fll.5p - fll.9l] . respectively, with 


initial data and (wg^^ Ug^^ tJ g\ g^^g). Denote by 

{6u,6v,6T,,6q,) = («('),n(i),TW,gW) - 


and 


(( 5 m 0 , SVq, STefi, Sqefl) 


[U, 




’ -‘■efly Hefli 


- [U 


( 2 ) 



Then, similar to fl3.8p . we have 


j^USu,h,,ST,M\\l + ^{\\Viu\\l + ||V*.||=) 

+ |5Tgp + \6qe\^ + |VTf)||5M||5Tg| + \Vqi^^\\6u\\5qe\]dxdy, (3.30) 


for all t G (0,7^. All the integrals on the right-hand side of the above inequality, 
except the last two terms, can be dealt with in the way as before in fl3.9p . while for 
the last two terms, we estimate them by the Holder, Sobolev embedding and Young 
inequalities as follows 


C [ {\VTp^\\6u\\6T,\ + \Vqi^^\\6u\\6qe\)dxdy 
7r2 

<C||Vri2)||^||5M||^||5Tg||2-FC||Vgf |U||hM||^||5ge||2 

m — 2 m — 2 
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+ c(||vrf>||^ + l|v«?»||^)(||ir.||2 + iiifciil), 

for alH e (0, 7^. Therefore, we deduce from fl3.30p that 

^11 ((5ti, (5 k, (5T„ (5(j„) III + i II ((Sii, (5ti) IlL 
< C(l + ||(KW,KW)||J + ||(Vt.W,VKW,VK«)||l)||(,5K,,5K,5T„,5,j.)||| 
+C(||VrP||,|, + ||V9P|||,)(||,5r,||| + 11,59,111), 
for all t E (0, 7~]. Applying the Gronwall inequality to the above inequality yields 

sup \\{6u,6v,6Te,6qe){s)\\l + ^ [ \\{6u,6v)\\Hids 

0 <s<t o Jo 

< gC'/j(l+||(A2)y2))||4 + ||(V«(2),Vd2),Vdl))||2 + ||(VrF,Vg?))||^)ds 

x\\{6uo,6vo,6Tefl,6qe^o)\\l, 

for all t G (0,7”]. Recalling the regularities in (i) and (iii), the above inequality 
implies the continuous dependence of the strong solution on the initial data on [ 0 , 7”], 
for any arbitrary T G ( 0 , oo). This completes the proof. □ 

4. Global existence and uniqueness of the limiting system 

In this section, we prove the global existence and uniqueness of strong solutions to 
the Gauchy problem of the limiting system fll.12p - fll.18p : 

Proof of Theorem I J.M (i) The global existence and regularities. By Theo¬ 
rem [m for any positive e, there is a unique global strong solution v^, qee) to 
system fll.5p - fll.9l) . with initial data (mq) 1^05 ^e ,05 7e, 0)5 such that 

sup ('MWI + ||(u^,u„T,.,fc)(()iiy') + r(t^Ail + ||(„^,u,)iiy')<i( 

0<t<T \ £ J Jo \ ^ ) 

+ [ {\\{dtUe,dtVe,dtTee,dtqee)\\l + \\'^Ue\\oo)dt<C, 

Jo 

for any positive hnite time T, where G is a constant depending only on a, Q, T and 
initial norms || (mq) P 05 ^e ,05 7 e,o)||Hi 5 and in particular, is independent of £. Moreover, 
if in addition that (VTe ^05 Vge,o) ^ for some m G (2,cx)), then we have 

further that 

sup \\{VTee,Vqee){t)\\l, < C {a,Q,T,m, ||(Mo,Uo,Te,o,ge,o)||Hi, || (VTe,o, Vge,o) ||m) , 
0 <t<T 

for any positive hnite time T, and, again, the estimate is independent of e. 
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Thanks to the above e-independent estimates, there is a subsequence, still denoted 
by {ue,Ve,Tee,qee), and {u,v,Te,qe), such that 

{Ue, Ve) ^ (m, n), in L^(0, T; 

{dtUe,dtVe) {dtu,dtv), in L^(0, T; T^(M^)), 

(Te„ gee) ^ {Te, qe), m L°°(0, T; 

(dtT,,, dtqee) {dtTe, Otqe), in L2(0, T; 

g+ ^ 0, in L-(0, T; p ^^(0, T; H\R^)), 

for any positive hnite time T, where ^ and ^ are the weak and weak-* convergences, 
respectively. The last convergence in the above implies that 

g^ = 0, or equivalently qe < 0, a.e. in x (0, T). 

Moreover, by the Aubin-Lions lemma, and using the Cantor diagonal argument, we 
have a subsequence, still denoted by {ue,Vs,Tee,qee), such that 

{ue,Ve) ^ Kn), in C{[0,r];L\BR))nL\0,T;H\BR)), 

(Tee, qee) ^ (Te, qe), in C([0, U L\Br)), 

for any positive hnite time T, and disc Br C of arbitrary radius i? > 0. 

Thanks to the previous convergences, one can take the limit £ —>■ 0+ in the equa¬ 
tions fll.5p - fll.8p for {ue,Ve,Tee,qee) to deduce that {u,v,Te,qe) satishes equations 
fll.5p - fll.8p . a.e. in x (0, cx)), since R in the previous strong convergences is arbi¬ 
trary; and moreover, by the lower semi-continuity of the norms, the a priori estimates 
stated in Theorem 11.21 hold. In order to complete the proof of existence, we still need 
to prove that qe satishes inequalities fll.16p - fll.18p . Inequality fll.lZp has already been 
verihed before. While for fll.l6p . note that equation fll.Qp for qee implies that 

dtqee + Ue ' Vqee + (Q + «)V • Ue < 0, a.e. in X (0, cxo), 

from which, recalling the previous convergences, one can take the limit e —)■ O’*" to 
see that 

dtqe + u ■ Vge + (0 + • n < 0, a.e. in x (0, oo), 

which is fOHU . 

It remains to verify fll.l8p . To this end, let’s dehne the set 

0~ = {(x, t)|ge(x, t) < 0, X G t G (0, oo)}, 
and for any positive integers j, fc, /, we dehne 

qe{x,t) < -j,x e Bk,t e (0,/)| , 

where B^ C is a disc of radius fc, and j, k,l eN. Noticing that 

o- = U“ ur=i 
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to prove that (11.181) holds a.e. on 0~, it suffices to show that it holds a.e. on 
for any positive integers j, fc, 1. Now, let’s £x the positive integers j, k, 1. Recalling 
that gee —t Qe in C{[0,T]; for any positive time T and positive radius R, it 

is straightforward that gee —t ge in Therefore, there is a subsequence, still 

denoted by gee, such that gee —)■ ge, a.e. on By the Egoroff theorem, for any 

positive number p > 0, there is a subset of with |E^| < r], such that 

gee -)■ ge, Uniformly on \ Er,. 

Recalling the dehnition of this implies that for sufficiently small positive e, it 
holds that ^ ^ 

qee<qe + ^< < 0, ou O"; \ E^. 

As a result, by equation fll.Op for gee, we have, for any sufficiently small positive e, 
that 

:= dtqee + Ue ' Vgee + (Q + a)V ' Ue = 0, a.e. on \ E^. 

Noticing that 

Qe dtqe + U ■ Vge + {Q + a)v -v G, in L^(0, T; 

for any positive hnite time T, which in particular implies Ge G, in L‘^{Ojki \ Erj). 
Since Ge = 0, a.e. on Ojki \ we have G = 0, a.e. on Ojki \ that is 

dtqe + u ■ Vge + (Q + a)V • n = 0, a.e. on Qjki \ E^,. 

By Lemma [T3l this implies that the above equation holds, a.e. on and further 
on 0~, in other words, (11.181) holds. 

Therefore, {u,v,Te,qe) is a global strong solution to system fll.12p - fll.18p . with 
initial data (mq, T'e,o, ge,o), satisfying the regularities stated in the theorem. 

(ii) The uniqueness. Let {u,v,Te,qe) and {u,v,Te,qe) be two strong solutions 
to system fll.l2l) - fll.l8|] . with the same initial data (mq, Uq, Te^o, ?e,o)- Define the new 
functions 

{6u, 6v, STe, Sqe) = (u, V, Te, ge) - (m, V, Te, ge). 

Then, one can easily check that {6u,6v,6Te,6qe) satisfies equations fl3.2p - fl3.5p . and 
the same argument as that for fl3.7l) yields 

<C [ [(|Vh| + |Vh| + |Vn| + |np + + l^np) 

+ \6Te\^ + \6qe\^ + \Vfe\\6u\\5Te\]dxdy. (4.1) 

We need to estimate dqe- To this end, we first derive the equation for 6qe- We 
divide the domain D := x (0 , cxd) as follows 


D = U Da U Ds U D4, 
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where 


= Ue < 0} n {qe < 0}, ^2 = {qe < 0} fl {q^ = 0}, 

^3 = {qe = 0} n {qe < 0}, ri4 = {qe = 0} n {qe = 0}. 

On the set Oi, qe and qe satishes, respectively 

dtqe + u ■ Vge + (Q + a) V • n = 0, 

dtqe + u ■ Vqe + {Q + a)V • h = 0. 

Subtracting the above two equations yields 

dtSqe + u ■ V5qe + 5u ■ Vqe + (Q + = 0, on Oi. (4.2) 

On the set 02, qe satishes 

dtqe + u ■ Vqe + (Q + a)V • n = 0, 

while for qe, since ge = 0 on O2, one has {dtqe, Vqe) = 0, a.e. on O2, and thus dtqe + 
-u • Vge = 0, a.e. on O2. Here, we have used the well-known fact that the derivatives 
of a function / G hh[oc(0) vanish, a.e. on any level set {{x,y,t) G Q:\f {x,y,t) = c}, 
see, e.g., 0 or page 297 of 0]. We will used, without any further mentions, this 
fact several times in the proof of this part. Therefore, one has 

dtSqe + u ■ V6qe + 5u ■ Vqe + (Q + a)V • n = 0, a.e. on O2. (4.3) 

Similar to (14.3p . on the domain O3, one has 

dt5qe -l- u ■ V6qe + Su ■ Vqe — {Q + a)V -7 = 0, a.e. on O3. (4.4) 

Finally, since ge = ge = 0, on O4, one has 

dtSqe + u ■ V6qe + Su ■ Vqe = 0, a.e. on ^4. 

Thanks to the last equation, as well as (I4.2p - (l4.4p . we obtain the equation for Sqe as 
dtSqe u ■ VSqe + Su ■ Vqe = -{Q + a)[V ■ Svxut + V • vxn 2 - V • vxqs] 

= -{Q + a)[V -Sv-V ■ Svxn^ + V • vxn 2 - V • vxus], (4.5) 

a.e. on X (0, cxd). Moreover, equation (14.51) holds in cxd); L^(R^)). 

Multiplying equation (14.51) by Sqe, and integrating over R^, then it follows from 
integration by parts that 
1 d 


2 dt 


ll'^?e||2 = ~ / [Su ■ VqeSqe + {Q + a)V ■ Sv]Sqedxdy 


-{Q +a) (V • vxn2 - V • vxn3){qe - qe)dxdy 


A 


\VSv\^dxdy + C / {\Sqe\^ + \Vqe\\Su\\Sqe\)dxdy 


{Q + a) (V • vxn2 - V • vxn3){qe - qe)dxdy. 


( 4 . 6 ) 
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Recalling that ge = 0 on ^ 2 , we have dtqe + u-Vqe = 0, a.e. on ^ 2 , and thus it follows 
from fll.lhp for (h, v, T^, qe) that V • h < 0, a.e. on ^ 2 - Similarly, one has V • n < 0, 
a.e. on ^ 3 . Thanks to these facts, we deduce 


V ■ vxn^iqe - ge) = V • vxn2qe > 0, 

-V • vxn^fle - ge) = V • vxn^qe > 0. 

Therefore, it follows from (14. 6 p that 

4 ll<^ge ||2 < aI|V 5 t ||2 + C f (l<^geP + \Vqe\\6u\\6qe\)dxdy. 

at / j ^2 

Summing the above inequality with fl4.1l) yields 

+1(IIVHI2 + Iivit>|l2) 

< c[ [(|Vm| + |Vt| + |Vt| + |Tp + |h|2)(|5M|2 + |5t;|2) 

iR 2 

+ |5Te|2 + |<5ger + |VTe||<5n||5Te| + \Vqe\\6u\\6qe\]dxdy, (4.7) 

which is exactly the same as inequality f|3.8p . from which, by the same argument as 
that in the proof of the uniqueness part of Proposition 13.11 one obtains 

\\{6u,6v,6Te,6qe)\\l = 0 . 


This proves the uniqueness. 

(iii) Continuous dependence. Let {vtS'^\v^'''\Te\q^e'^) be the unique solutions 
to system fll.l2l) - fll.l8l) . with initial data {uq\vq\t^^q, z = 1 , 2 . Suppose, in 
addition that (VT^ q, Vgg*o) G L™(M^), for some m e (2, 00 ). Then, recalling what we 

have proven in (i), {u^^\v^^\Te^\qe'^) has the additional regularity that {Te\qe'^) G 
T°°(0, T; L™'(R^)), for any positive time T. 

Denote by 


{Su,Sv,ST,,Sq,) = (u<‘>,i,<‘>,ri‘>,gW) - 


and 


{ 6 Uo, 6 Vq, STefi, Sqefl) 


[u, 




7 ^( 1 ) „(!)' 

' -^e ,0 5 ge ,0 2 


- [U 


( 2 ) 



Then, similar to fl4.7p . we have 


j^\\{Su,Sv,ST„SqMl + l(l|Vfe||^ + llVfol© 

< cf [(|Vu(2)| + + |Vud)| + |z;d)|2 + |^;(2)|2)Q^^|2^ 

+ |5Tg|2 + |5gg|2 + \VTP^\\5u\\ 6T,\ + \Vqi^^\\6u\\5q,\]dxdy, 
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which is exactly of the same form as (I3.30p . Therefore, by the same argument as that 
in the proof of the continuous dependence part of (iii) of Theorem 11.11 we obtain 

sup \\{6u,6v,6Te,6qe){s)\\l + ^ [ \\{6u,6v)\\Hids 

0<s<t o Jq 

< gC/o*(l+||(«(2),^(2))||| + ||(V«(2),Vd2),Vdi))||2 + ||(VrP,Vg?))||^)<is 

x\\{6uo,5vo,STefl,6qe,o)\\l- 

Recalling the regularities of {u^^\v^^\Te\qe^), i = 1,2, the above inequality implies 
the continuous dependence of strong solutions on the initial data. This completes 
the proof of Theorem 11.21 □ 


5. Strong convergence oe the relaxation limit 


In this section, we prove the strong convergence of the relaxation limit, as £ —)■ O’*", 
of system fll.5l) - fll.9p to the limiting system fll.12p - fll.18p : 

Proof of Theorem \1.31 Dehne the difference function {6ue, Svs, ST^e, 6qee) as 

{5u^, 6v^, 6Tee, Sqee) = («£, V^, T^e, qee) “ («, V, Tg, qf). 

Taking the subtraction between equations fll.SD - fll.Sp . for (ue, Ue, gee), and equa¬ 
tions fll.i 2 p - fll.i 5 p . for {u,v,Te,qe), one can easily check that 


dt5u^ + {6us ■ V)Sus + {due ■ V)u + {u ■ V)6us — Adu^ 

+ V5p^ + V ■ {5v^ 0 5v^ + 5v^ ®v + v® 5v^) = 0, 
V ■ 5us = 0, 

dtdVe + {5Ue ■ V)5Ve + {5Ue ■ V)n -|- {u ■ V)5Ve — A5Ve + {5Ve ■ V)5Ue 


-f- {5Ve ■ V)m -I- {v ■ V)5Ue = 


1 -f- o 


V(5Te, - 5q,,), 


dtSTee + 5 Ue ' V5Tee + ' VTe U ■ VST^e “ (1 “ Q)V • = 0, 


(5.1) 

(5.2) 

(5.3) 

(5.4) 


where fl5.ip - fl5.4p hold a.e. on x (0, cx)) and in LfQ^([0, cx)); L^(R^)). 

Multiplying equations fl5.ip . fl5.3p and f|5.4p by Su^, Sv^ and ST^e, respectively, 
summing the resultants, integrating over R^, and noticing that 

/ [V ■ 0 6vi;) ■ 6ui; + {6v^ ■ V)5u^ ■ 6vi;]dxdy = 0, 

Jr 2 

it follows from integration by parts that 


1 d 

2 dt 


\\{6u„6v„6T,,)\\l + \\{V6u„V6v,)\\l 


[{Sus ■ V)m -|- V ■ {6ve ®v + v® • du^dxdy 
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[{6us ■ V)t + {6vs ■ V)m + {v ■ V)(5m£] • dv^dxdy 


1- a 


(V • 6vs){6Tee - 5qee)dxdy 


- / [6ue ■ VTe - (1 - Q)V • 6ve]6Teedxdy, 


from which, by the Young inequality, we deduce 


]^d\\(Su,,6v„ST„)\';i+\\(V6u„VSv,)\\l 

< [ [(|Vm||(5m£| + 2|Vu||(5u£| + 2|u||V(5u£|)|(5m£| + (iVull^Mel 

2R2 

+ \'\/u\\6ve\ + \v\\'\/6ue\)\6ve\]dxdy + - [ IVdVelilSTeel + \Sqee\)dxdy 

1 +a J^2 


+ [ [(1 - Q)\V6v,\\6Tee\ + \VTe\\5u,\\6T,,\]dxdy 
Vr2 

<11 {\V6u,\^ + \V6v,\^)dxdy + C [ [(|Vu| + |Vu| + |up) 

2 Jr2 Jr2 

X{\Su,\^ + + \6Tee\^ + \6qes\^ + \VTe\\6Us\m,\]dxdy. 


Therefore, we obtain 

j^\\{Su,,Sv,,SnM + \\i^S^s,VdvM 

< C [ [(|Vu| + |Vu| + |un(|5u,|2 + |5^;,|2) 

7r2 

+ |5Te,|2 + \dqee\^ + | VT^ 11115Te, |]da:di/. (5.5) 


We still need to estimate H^q'eelli- To this end, we first derive the equation for 6qee- 
On the set {{x,y,t) G x {0,oo)\qe{x,y,t) < 0}, gee and ge satisfy equations fll.91) 
and fll.181) . respectively, and thus 5gee satisfies 


dtdqee + 6ue ■ V6qee + Sus ■ Vge + u ■ VSqee + (Q + a) V • 6Vs 


1 H“ 


-Qee^ 


a.e. on {{x,y,t) G x (0, cxD)|ge(a;, g, f) < 0}. On the set O := {(x,t) G x 
(0, cxd) |ge(a;, t) = 0}, recalling, again, the well-known fact that the derivatives of 
a function / G 1T[q(!(R^ x (0, cxd) vanish, a.e. on any level set {{x,y,t) G x 
{0,oo)\f{x,y,t) = c}, we have dtqe + u ■ Vge = 0, a.e. on O, and gee satisfies fll.Qp . 
Consequently, 6qee satisfies 


dtSqee + dUe ■ V6qee + 5Us ■ Vge “h U ■ V6qee + {Q + a)V • Te 


1 -l- tt I 
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a.e. on O. Combing the above two equations, one can see that 5qee satishes 

dt5qee + SUe ■ ^Sq^e + SUe ' V^e + U ■ V6qee 
— 1 “h — 

+ (0 + a)V ■ 6 vc = -gj - (<3 + a)V ■ vxo{x, y, t), 


(5.6) 


a.e. on x (0, cx)) and in cx)); L^(R^)). 

Multiplying equation fl5.6p by Sqee, and integrating over R^, then it follows from 
integration by parts that 


1 d 


1 -|- Q) 


2 A £ 


qteSqeedxdy 


= - [6ue ■Vqe +{Q + ■ 6ve + V ■ vxo{x,y,t))]6qeedxdy, 

a.e. t G (0, CX)). Recalling that qe < 0, we have 

f qiSq„dxdy = f gj(g„ - q,)dxdy > f qU^dxdy = IlgJU^. 


(5.7) 


(5.8) 


Note that dtqe + u ■ Vge = 0, a.e. on C>, it follows from fll.161) that V • n < 0, a.e. on 
C>, and thus 

- V ■ vxo{.x, y, t)5qee = -V • vxo{.x, y, t)qee < -V • vxo{.x, y, t)q+. 

Thanks to the above inequality, it follows from fl5.7l) . fIS.Sp and the Young inequality 
that 


< 


[Sus ■ Vge + (Q + «)V • 6vs]Sqeedxdy — {Q + a) V • vxo{x, y, t)q^^dxdy 


< / \Vqe\\Sue\\6qee\dxdy + -\\V6ve\\l + {Q + afWSq, 

./n2 4 


e £\\2 


1 -|- 

and thus 
d 


hts\\l + 


dt e 


{Q + Qc)^ 
2(1 + a) 

2 , 1 “f" 




Ike^£ll2 <^l|V5us||2 + 2(Q + aY\\6q^,\\l + ^ g||Vu||^ 


(Q + aY 


1 “h Q; 


\Vqe\\due\\5qee\dxdy. 


(5.9) 


Summing fl5.5p with fl5.9p yields 
d 


1 -|- 


^^\\i6u,M:5T,,,6qeem + -||(V5u„ V(5u,)||^ + 


Ik, 


+ ||2 
e£ II2 
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<6* / [(|Vm| + |Vt| + \v\^){\Su,\^ + 

+ \VTe\\6u,\\6Tee\ + \Vqe\\6u,\\6qee\]dxdy + Ce\\Vv\\l 
from which, it follows from the Hdlder, Ladyzhenskay, Gagliardo-Nirenberg, ||(/7|| 2 m < 

m — 2 

m — 2 2 

^Il‘dll 2 ’" Young inequalities that 

±\\(5u,,5v,,5T,,,5q,M + ^\\{V5u,,V5v,)\\l + ^ 

<C (II (Vn, Vn)II 2 + ||n||2) || (5n„ Sv^) ||^ + G|| Sq,,) 

+ C'||(VTe,Vge)IUI|5M.||AH^||(5Te„<5ge.)||2 + C'£||Vn||2 

m—2 

<C'(||(Vn,Vn)||2+||n||2||Vn||2)||(<5n.,<5n,)||2||(V<5n„V5n,)||2 + C'£||Vn||2 

m — 2 2 

+ C'||(<5Te„<5ges)||2 + C'||(VTe,Vge)|U||5we|l2'" IIV<5n,|| 2 "* ||(<5T,„ 5^^.)||2 

<^||(V5n,,V<5n,)||2 + G(||(Vn,Vn)||2 + ||n||2||Vn||2) \\{Su,,Sv,)\\l + Ce\\Vv\\l 

m m — 2 m 

+ C||(iT„,ife)||l + C||(VT„V9.)ll™-|l'5«,ll2“ll('5T„,ife)||r-‘ 
<i||(Vi«.,Vi't>.)ll2 + c(l + ||(V«,V«)ll2+ll'o|l2l|Vw|l^ 

+ ||(VTe,Vge)llm-^)||(5w.,5t^.,5Te,,<5ges)|l2 + C'£||Vn||2. (5.10) 

Therefore, we have 

( m \ 

1 + ||(Vn, Vn)||2 + ||n||2||Vn||2 + ||(VT„ Vge)||^-^) 

X ||(<5n„5t;„<5Te„5ge.)||2 + C'c||Vn||2. 

Applying the Gronwall inequality to the above inequality and recalling the regu¬ 
larities of {u,v,Te,qe) yield 

sup \\{5u,,5v,,5T,,,5q,,){t)\\l+ (\\(V5u,,V5vM+^^^) dt < Ce, 

0<t<T Jo \ ^ / 

for a positive constant C depending only on a,Q,T,m, || (wo; T), q, ?e,o)||iyi and 
||(VTe^ 0 ; Vge,o)||m- This proves the desired estimate in the theorem, while the strong 
convergences are direct consequences of this estimate. □ 

6. Appendix 

In this appendix, we state and prove several parabolic estimates, which have been 
used in the previous sections. 
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Lemma 6.1. Given a time T G (0,cx)), and a function g G L"(0, T; with 

1 < a, (3 < oo. Let U be the unique solution to 

dtU-AU = g, znM2x(0,r), 

U\t=o = 0, in 

Then, we have the estimate 

ll^ff^llL“(0,r;L'5(R^)) + ll^^ll-L“(0,r;L^*(R2)) < Ca^/3\\g\\L°‘(0,T-,L^(^^))^ 

where is a positive constant depending only on a, ft, and in particular is inde¬ 
pendent ofT and g. 

Proof. Introducing the scaled functions Ur and gr as 

Ur{x,t) = U{VTx,Tt), gr{x,t) = g{VTx,Tt), a; G R^, t G (0,1), 

then one can easily verify that Ur and gr satisfy 

dtUr-AUr = Tgr, inR2x(0,l), 

U\t=o = 0, in R^. 

Applying the maximal regularity theory for parabolic equations to the above system 
(see, e.g., [l3[, jl^ and jlsj), one has 

WdtUrh “(0,l;L^(R2)) + ll^tVllL“(0,l;L/3(K2)) < Ca,j3T\\gr\\L°‘{0,T;Ll^{R^))- 

From which, and after observing that, 

||^thV||L“(0,l;L/3(K2)) = “ ^||c?t[/|| L“(0,r;L^(K2)), 

ll^fVllL“(0,l;L/3(]R2)) = “ ^||At/||ia(0,r;L'3(K2)), 

ll5'rllL“(0,l;L/’(R2)) = “ ^||fi'||L“(0,r;L^(R2))5 

one obtains the conclusion. □ 

Lemma 6.2. Given a time T G (0, cxd), and let f and g be two functions, such that 
f G L^(R^ X (0,T)) and g G L"‘(R^ x (0,T). Let v be the unique solution to 

dtv — Av = f + Vg, in R^ x (0, T), 
v\t=o = Vo ^ m R^. 

Then we have the following estimate 

2+^y Il/ll2^y + J ) 

where C is an absolute constant, and in particular is independent ofT,vo, f and g. 
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Proof. Decompose v as v = v + v, where v and v are the unique solutions to systems 

f dtv — Av = /, in x (0, T), 

1 v\t=o = Vo e in R^, 


and 

jdtV-Av = Vg, inR2x(0,T), . ^ 

\ v\t=o = 0, in R2, 

respectively. The standard energy approach (multiplying the equation for v by —Av, 
integrating over R^, integration by parts, using the Young, and integrating with 
respect to t over (0, T)) to the system for v leads to 


sup ||VT(t)|| 2 + [ ||Ar;|| 2 dt < IIVT 0 II 2 + 
o<t<r Jo 



WfWldt. 


Dehning U to be the unique solution to the system 

{dtU-AU = g, inR2x(0,r), 

\ U\t=o = 0, in RT 

Then Vf/ satishes the same system as that for v, and therefore, by the uniqueness 
of the solutions to system fl6.1l) . we have v = VU . Thanks to this fact, and applying 
Lemma 16.11 it follows from the elliptic estimates that 


ll^'^^l|L4(o,r;L4(R2)) — II V^f/||2,4(0^7-;2,4('R2)) 

< C|| Af/||i4(o,r;L-‘(lR2)) < C||5f||2,4(0,r;L4(lR2)), 

for an absolute positive constant C. 

Combining the estimates for v and v, we deduce from the Ladyzhenskaya inequality 
that 

r w^vwtdt <c r\\vv\\idt+c r\\vv\\tdt 

Jo Jo Jo 

<c (||vi.„|i5 + ii/iipi) + iisiipi 

for an absolute positive constant C. This completes the proof. □ 

Lemma 6.3. Given a time T G (0, 00 ) and a number m G (2, 00 ). Let f G 
L^(0, T; L"*(R^)), and v be the unique solution to 

f dtv — Av = /, in R^ x (0, T), 

\ v\t=o = VoeH\RP). 
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Then, we have the following estimate 

rT 


||Au||„irfi < C„t{l + VT) 


IIVm 


0||2 


r-r \ 5'' 

|2 


11/11™^^ 


where Cm is a positive constant depending only on m, and in particular is independent 
ofT,f and vq. 

Proof. Decompose v as v = v + v, where v and v are the unique solutions to systems 

dtv — Av = /, in x (0, T), 
v\t=o = 0, 

dfV — Av = 0, in R^ x (0, T), 
v\t=o = VoeH\R^), 


and 


( 6 . 2 ) 


respectively. 

By Lemma [6. II and using the Holder inequality, for v, we have the estimate 


r-r 


\\Av\\mdt < T 2 


rT 


\\Av\\idt] <CmT-^ 


r-r 


wnidt 


To estimate v, we multiplying equation (I6.2p by tA'^v — Av, integrating the resultant 
over R^, then it follows from integration by parts that 

~(l|v«||^ + WViMWl) + 1|| A«||l + ||yiVA«||^ = 0, 

from which, integrating with respect to t yields 


rT 


sup (||Vh(t )||2 + \\ViAv{t)\\l) + / (||Ah||2 + \\ViVAv\\l)dt < ||Vno| 


o<i<r 


Thanks to this estimate, by the Gagliardo-Nirenberg, 
Holder inequalities, we deduce 


1-^ 


< C'llv?ll 2 "*l|V (^||2 ™,and 


rT 


rT 


\\Av\\mdt<C I ||Ah||^||VA {)||2 '"dt 


r-r 


=C / II Ahll^ ||\/fVA'0||2 "^t 2^^ ^^dt 


r-r 


IIAhll^df 


<G 


<Cy/rnT II Vuolb- 


r-r 


WViV Av\\ldt 


m — 2 
2m 


rT 


t ^'^dt 
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Combining the estimates for v and h, we then deduce from the Young inequality 
(recalling m > 2) that 

“T "T 

f \\Av\\mdt< [ {\\Av\\jn + \\Av\\m)dt 

Jo J 

T' 

>0 


'0 




dt] + Ci/mT™ IIVU0II2 


+ V^) 

proving the conclusion. 


IIVU 0 II 2 + / Wfltdi 


□ 
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